Baryon Acoustic Oscillations in 2D II: Redshift-space halo clustering in N-body 

simulations 
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We measure the halo power spectrum in redshift space from cosmological N-body simulations, 
and test the analytical models of redshift distortions particularly focusing on the scales of baryon 
acoustic oscillations (BAOs). Remarkably, the measured halo power spectrum in redshift space 
exhibits a large-scale enhancement in amplitude relative to the real-space clustering, and the ef- 
fect becomes significant for the massive or highly biased halo samples. These findings cannot be 
simply explained by the so-called streaming model frequently used in the literature. By contrast, 
a physically-motivated perturbation theory model developed in the previous paper reproduces the 
halo power spectrum very well, and the model combining a simple linear scale-dependent bias can 
accurately characterize the clustering anisotropies of halos in two dimensions, i.e., line-of-sight and 
its perpendicular directions. The results highlight the significance of non-linear coupling between 
density and velocity fields associated with two competing effects of redshift distortions, i.e., Kaiser 
and Finger-of-God effects, and a proper account of this effect would be important in accurately 
characterizing the BAOs in two dimensions. 



I. INTRODUCTION 

Redshift distortions caused by the systematic effect of 
peculiar velocity of galaxies induce the anisotropies in the 
galaxy clustering patterns (e.g., 0, Q)- These are now 
recognized as a powerful tool to test theory of gravity 
on cosmological scales with a great interest (e.g., 0-0] )■ 
On large scales, the magnitude of redshift distortions is 
simply characterized by the growth-rate parameter /, de- 
fined as / = dhiD + /dliia, where the quantities D + and 
a respectively denote the linear growth factor and the 
scale factor of the Universe [H, H Q. Since the modi- 
fication of gravity from the general relativity generally 
alters not only the cosmic expansion but also the struc- 
ture formation, the measurement of growth-rate parame- 
ter provides a complementary way to distinguish between 
various models of gravity. 

Notice that distortions of the galaxy clustering pat- 
tern also arise from apparent mismatch of the under- 
lying cosmology when we convert the redshift and an- 
gular position for each galaxy to the comoving radial 
and transverse distances. This is known as the Alcock- 
Paczynski effect [ljj, and with the baryon acoustic os- 
cillations (BAOs) as a standard ruler, it can be utilized 
for a measurement of both the Hubble parameter H(z) 
and angular diameter distance Da(z) of distant galax- 
ies at redshift z (e.g., fllTfl5l |. see [l6l - l24| for observa- 
tional status). Hence, a measurement of the anisotropic 
galaxy clustering serves as a dual cosmological probe 
to simultaneously constrain both the cosmic expansion 
and structure growth (e.g., [25]). Planned and ongoing 
galaxy redshift surveys such as the Baryon Oscillation 
Spectroscopic Survey (BOSS) [H, Hobby-Eberjy Tele- 
scope Dark Energy Experiment (HETDEX) [27], Sub- 
aru Measurement of Imaging and Redshift equipped with 



Prime Focus Spectrograph (SuMIRe-PFS) H, and EU- 
CLID/ JDEM [29, [30] aim at precisely measuring the 
anisotropic power spectrum and/or two-point correlation 
function in redshift space, from which we can simultane- 
ously determine Da, H and / in a very accurate way. 

To get tight and robust cosmological constraints, an 
accurate theoretical template for the anisotropic power 
spectrum is highly demanding, taking a proper account 
of various systematic effects. In particular, in redshift 
space, in addition to the non-linear gravitational clus- 
tering, the clustering statistics generally suffer from two 
competing effects of redshift distortions, i.e., enhance- 
ment and suppression of the clustering amplitude, re- 
ferred to as the Kaiser and Finger-of-God effects, respec- 
tively [H [3~l| . While the Kaiser effect comes from the co- 
herent motion of matter (or galaxies) and the magnitude 
of this effect is simply described by the growth-rate pa- 
rameter /, the Finger-of-God effect is mainly attributed 
to the virialized random motion of the mass residing at 
a halo, and the significance of this is rather sensitive to 
the properties of the small-scale clustering. In the weakly 
non-linear regime, a tight correlation between the veloc- 
ity and the density fields still remains, and a mixture of 
Kaiser and Finger-of-God effects is expected to be sig- 
nificant. In this sense, a careful treatment is needed to 
model the anisotropic power spectrum accurately, other- 
wise one might derive a biased estimate of the growth- 
rate parameter as shown by, e.g., [32j . 

In a previous paper [33[ , based on the analytical treat- 
ment with perturbation theory, we have presented an 
improved prescription for redshift-space power spectrum 
relevant for the scales of BAOs. The model properly 
takes account of both the non-linear gravitational clus- 
tering and redshift distortions. Contrary to the so-called 
streaming model (e.g., @, HH), which has been phe- 
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nomenologically introduced but frequently used in the 
literature, the new model includes the corrections coming 
from the non-linear coupling between velocity and den- 
sity fields, which gives rise to a slight increase in the am- 
plitude of the redshift-space power spectrum. Combining 
the improved treatment of perturbation theory to accu- 
rately predict non-linear clustering [35l [36j , the model re- 
produces the monopole and quadrupole moments of the 
matter power spectrum measured from N-body simula- 
tions quite well, and a percent level precision is achieved 
over the scales of interest for BAOs 33]. However, the 
comparison with N-body simulations has been only made 
with lower multipoles of the matter power spectrum, i.e., 
monopole and quadrupole spectra. 

In this paper, as a natural extension of the previous 
study, we further test the analytical model of redshift 
distortions. Using a large set of N-body simulations, 
we measure the redshift-space power spectrum in two 
dimensions, characterized as a function of line-of-sight 
wavenumber k\\ and its normal one k±. In particular, we 
examine the halo power spectrum in detail, and investi- 
gate the extent to which the perturbation theory descrip- 
tion combining a simple halo bias scheme can describe the 
halo clustering properties in redshift space. We find that 
the measured halo power spectrum in redshift space at 
relatively large scales k < 0.2 /iMpc^ 1 exhibits a signifi- 
cant enhancement in amplitude relative to the real-space 
clustering, which cannot be explained by the stream- 
ing model. On the other hand, the new model includ- 
ing the corrections can describe the enhancement fairly 
well. The results indicate that the halo bias can illumi- 
nate the non-linear coupling between density and veloc- 
ity fields, and a proper account of the non-linear velocity- 
density coupling seems important in accurately modeling 
the galaxy/halo power spectrum in redshift space. 

This paper is organized as follows. We begin with a 
brief review of the models based on perturbation theory 
in redshift space in Sec. [TTJ Sec. IIIII describes the setup 
of N-body simulations and discusses the power spectrum 
analysis. The results for a detailed comparison between 
models and simulations are presented in Sec. HVl particu- 
larly focusing on the halo and matter power spectrum in 
two dimensions. We discuss about the impact of the non- 
linear correction terms on some typical on-going/planned 
galaxy redshift surveys in Sec. El Finally, Sec. [VI] is de- 
voted to the summary of this paper. 



II. THE MODELS 

In principle, the clustering statistics in redshift space 
can be mapped from those in real space through the coor- 
dinate transformation between the positions in real and 
redshift spaces, r and s 66]: 



where the quantity H is the Hubble parameter at red- 
shift z, the unit vector z indicates the line-of-sight di- 
rection, and the quantity v z represents the line-of-sight 
component of the peculiar velocity field. Although the 
expression ([T]) is very simple, the statistical relationship 
between the real- and redshift-space quantities is rather 
complicated and is actually difficult to treat without ap- 
proximations. One reason for the difficulty comes from 
the anisotropics induced by the velocity field which is 
very sensitive to the small-scale structure. The other im- 
portant aspect of the redshift-space clustering is the cou- 
pling between velocity and density fields, which produces 
an apparent structure growth by the coherent motion. 

Qualitatively, the clustering properties in redshift 
space are explained by the Kaiser and Finger-of-God ef- 
fects. The so-called streaming model is a phenomeno- 
logical model that accounts for these two effects sepa- 
rately @, [34]]. Here, among various streaming models, we 
consider the following form of the redshift-space power 
spectrum (e.g., [37l [3§]): 



P(k,(i) = D { (knfa v ) 

X [Ps 5 (k)+2ftl 2 P S g(k) 



f^P ee (k)], (2) 



i 



H{z) 



u«(r)z, 



(1) 



where / is the growth rate parameter, the quantity a v is 
the one-dimensional velocity dispersion, and /1 is the di- 
rectional cosine of the angle between line-of-sight z and 
the Fourier mode k. In the above, while the function 
Df represents a damping function mimicking the Finger- 
of-God effect, the term in the bracket indicates an im- 
proved prescription of the Kaiser effect which takes ac- 
count of the nonlinear gravitational evolution [iirj . The 
spectra P$$, Pge, and P$e respectively denote the auto- 
power spectra of density and velocity divergence, and 
their cross-power spectrum, with the velocity divergence 
9 defined by 6 = -(I + z)/(Hf)\7v. The explicit func- 
tional form of the damping function will be specified 
later. 

Notice that the streaming models have been originally 
introduced and frequently used in the literature to ex- 
plain the observed power spectrum on small scales. On 
large scales of our interest, non-linearity of the gravita- 
tional clustering is rather mild, but a tight correlation 
between density and velocity fields still remains. In pre- 
vious paper [33], starting with the rigorous expression 
of redshift-space power spectrum, we partially applied 
a \ow-k expansion, and found that the model ([2]) misses 
some important terms. The missing terms naturally arise 
from the next-to-leading order corrections of the low-fc 
expansion, and they represent the nonlinear coupling be- 
tween velocity and density fields. The model is given by 

P(k, t i)=D i {k^fa v ) 

x P ss (k) + 2 fti 2 Pse(k) + f^Peeik) 

+A{k,Kf) + B(k,mfj\. (3) 
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The newly derived terms A and B are written as 

/,;/) = (*/,/) 

x {B CT (p, k p, -k) - S CT (p, k, -k - p)} , (4) 

B(k, (i; f) = (fc/i ff J ^ F(p)F(k - p) ; (5) 
F(p) = p|p5 fl (p) + /^P M (p)| ) 
where the function _B CT is the cross bispectra defined by 



fl(kO ^ 5(k 2 ) + / -^(k 2 ) j |<5(k 3 ) + / -^f #(k 3 ) 

(2^) 3 <5 D (k 1 + k 2 + k 3 ) B a (k x , k 2 , k 3 ). (6) 



Within the range of the validity of the \ow-k expansion, 
these corrections should be small, and we apply the stan- 
dard PT treatment to calculate the A and B terms. To 
compute the corrections, we use the expressions (A3) and 
(B4) presented in Ref. [33| , which are suited for a fast nu- 
merical calculation. 

The previous study [33| reveals that the corrections 
slightly enhance the power spectrum amplitude over the 
scales of BAOs, and the model $Ji§ excellently reproduces 
the monopole and quadrupole spectra obtained from N- 
body simulations. Note, however, that the aforemen- 
tioned models of redshift distortions are those for the 
matter distribution, and the previous study restricted 
the analysis to the case of the matter power spectrum. In 
this paper, adopting a simple linear bias relation without 
no velocity bias in real space, i.e., Sh = &<5 m , we further 
test the model ^ against the halo clustering in redshift 
space. Then, the model of halo power spectrum becomes 

Ph(k,n) = D{(k[if<j v ) 



x b 2 



P SS (k) + 2^ 2 P Se {k) + P 2 H 4 Pgg(k) 

+bA(k,n;P) + b 2 B(k,n;P) (7) 



with the quantity (3 defined by f3 = f/b. We use the 
improved PT treatment developed by Refs. [HI, HH to 
compute Pss, Pee and Pse- We will check the validity 
and accuracy of the model prescription with and without 
the corrections A and B in detail for the halo cluster- 
ing. In what follows, we specifically call Eq. ([7]) includ- 
ing the corrections the TNS model, and discriminate it 
from Eq. ([7]) neglecting the corrections as the streaming 
model. As for the damping function Df, we adopt the 
Gaussian and the Lorentzian forms often used in the lit- 
erature (e.g., |39l442| ). and compare between the results 
of these predictions: 



D { (x) = 



exp(— x 2 ), 
1 

{l + x 2 /2f 



(8) 



Since the damping function mainly alters the global 
shape of the power spectrum and our primary focus is 
to model the acoustic feature precisely, we will treat the 
velocity dispersion <r v in the damping function as a free 
parameter, and determine it by fitting the prediction to 
the N-body simulations. 

Finally, note that in the case of the TNS model, the 
expression ([Jj is valid only when the linear bias param- 
eter b is scale- independent. Nevertheless, in later analy- 
sis, we allow to incorporate the scale-dependence of the 
bias into the model (UJ), and compare it with the halo 
power spectrum in redshift space. Strictly speaking, the 
scale-dependence of the bias changes the structure of the 
integral kernel in Eqs. (@| and ([5]), and we cannot use 
the formulas for the corrections A and B presented in 
Ref. [33| . which have been derived for the matter power 
spectrum. As shown in Sec. IIV A[ however, the scale- 
dependence of the halo bias calibrated from the N-body 
simulations turned out to be very weak over the scales of 
our interest. Thus, our treatment of the model (J7J with 
scale-dependent bias would be validated, and does not 
change the final conclusion. 



III. ANALYSIS 

A. N-body simulations 

To assess the validity and accuracy of the analytic 
models against the redshift-space halo clustering, we run 
a set of N-body simulations and identified halos. All the 
simulations are performed with a publicly available tree 
PM code GADGET2 [43|. We adopt N = 1,280 3 particles 
in boxes with a side length of 1, 144.72 h~ x Mpc [i.e., the 
volume is 1.5 (/i _1 Gpc) 3 ], and set the softening length to 
5% of the inter-particle distance. We generate the initial 
conditions at z = 99 by second-order Lagrangian pertur- 
bation theory (e.g., [43. l45j|) starting from particles put 
on the regular lattices. The matter transfer function is 
calculated by a publicly available Boltzmann solver CAMB 
[4^ | , assuming the best- fit ACDM cosmological model de- 
termined by the five-year WMAP [ItJ • The output data 
of 15 independent realizations are stored at z = 0.35. 
Table U summarizes the settings of the simulations. 



B. Halo catalogs 

Given the data set of dark matter clustering, we iden- 
tify the halos by a Friends-of-Friends finder with linking 
length 0.2 times the mean inter-particle distance. For 
each realization, we construct nine halo catalogs with dif- 
ferent minimum and maximum masses. Allowing a slight 
overlap of the mass range, each mass bin is determined 
so that the signal-to-noise ratio of halo power spectrum 
measured at k = 0.0975/iMpc -1 becomes roughly com- 
parable among the nine halo catalogs [see Eq. ([12")) for 
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TABLE I: Summary of simulation parameters. 



# of runs 


# of particles 


box size 


softening 


mass/particle 




■^out 




fib / fim 


h 


n s 


C8 


15 


1,280 3 


1, 144.72ft _1 Mpc 


44.72/i _1 kpc 


5.54 x W 10 h~ 1 M @ 


99 


0.35 


0.279 


0.165 


0.701 


0.96 


0.817 



the definition of signal-to-noise ratio]. Detailed proper- 
ties of the nine halo catalogs, including the mass range 
and the number density of halos, are summarized in Ta- 
ble HU In what follows, we especially call the least and 
most massive halo catalogs (bin 1 and 9) as light and 
heavy, and discuss their clustering properties in detail. 
Note that the heavy catalog has roughly the same values 
of the bias parameter, number density, and survey vol- 
ume as observed in the luminous red galaxies (LRG) [48[ 
of the seventh data release of Sloan Digital Sky Survey 
(SDSS DR7) 



C. Power spectrum 

We measure the dark matter and halo power spectra 
from the N-body simulations using the standard method 
based on the fast Fourier transform. To be specific, we 
first assign the dark matter or halo particles onto a 1 , 024 3 
grid by cloud-in-cells (CIC) interpolation scheme [50| to 
obtain a density field on a lattice. We then transform 
it to Fourier space and divide it by the CIC kernel in 
order to eliminate the window effect. Finally, we take 
the average of the power over the modes in each k bin 
in real space, and in (k,fj,) bin in redshift space. The 
size of each Fourier bin is set to Afc = 0.005/iMpc _1 
and A/i = 0.05 for the wavenumber and the directional 
cosine, respectively. 

Our primary focus is to test the analytic models of 
redshift distortions against the halo clustering in N-body 
simulations. To do this, we need to incorporate the halo 
bias properties characterized by the function b(k) into the 
model prediction ([7]). While, in practice, the quantity 
b(k) should be parametrized by a simple function, and 
must be determined by fitting the model prediction to the 
observed power spectrum, in order to make a transpar- 
ent test, we directly measure the halo bias parameter b(k) 
from the N-body simulation, and use it to compute the 
analytic models of redshift-space power spectrum. Note 
that this treatment still includes several non-trivial as- 
pects. Indeed, the linear bias prescription seems rather 
simplistic, and it might not capture the real clustering 
nature of halos over the BAO scales. Further, in redshift 
space, the presence of correction terms in the TNS model 
implies that even the simple linear bias can lead to a non- 
trivial modulation of clustering amplitude, which might 
drastically alter the acoustic structure in redshift space. 

Let us discuss how to characterize and measure the 
halo bias properties from N-body simulations. According 
to the treatment by Ref. [5l|, we decompose the halo 



density contrast defined in real space into two pieces: 

<5 h (k) = b(k)5 ni (k) + c (k). (9) 

Here, the quantity e, which satisfies (<5 m e) = 0, can be re- 
garded as the residual noise contribution, and is related 
to the stochasticity of the halo sampling process. No- 
tice that the above decomposition does not rely on the 
assumption of linear bias. The other additional contribu- 
tions, which cannot be simply written as the term linearly 
proportional to 5 m , are all incorporated into the second 
term. Given the auto and cross power spectra of matter 
and halos measured in real space, P m , Ph and Pi im , the 
decomposition of the signal linearly proportional to the 
matter density field, and noise can be uniquely done, and 
the halo bias parameter can be measured in the following 
way: 

b(k) = P hm (fc)/P m (fc), (10) 
P £ (fc) = P h (k)-b 2 (k)P m (k), (11) 

where the quantity P e is the residual noise power spec- 
trum defined by (e(k)e(k')) = (27r) 3 <5 D (k + k')P e (k). In 
deriving the above relations, we have used the property 
(S m e) = 0. Note that the halo subsamples listed in Ta- 
ble fll] are constructed so as to take the same signal-to- 
noise ratio among the nine catalogs. The signal-to-noise 
ratio is defined by 

/S\ _ b\k)P m {k) 

UA~ ?^ ■ { } 

At k = 0.0975 /iMpc -1 , the signal-to-noise ratio roughly 
reaches (S/N)^ ~ 4 for all of our subsamples. 

In testing the analytic models of redshift-space power 
spectrum, we subtract the residual noise contribution 
from the measured halo spectrum, and then compare the 
results with analytic prediction in the (k, fi) space. Note 
that apart from the bias b(k), the models and © 
[or ([7])] involve a single parameter er v , which has to be 
determined by fitting the model power spectrum to the 
measured power spectrum. To derive the best-fit value 
of er v , we assume the Gaussianity, and use the statistical 
error of the redshift-space power spectrum given by [52j 

[AP(fc,/.)] 2 - -J— [b 2 (k)P m (k,ri + P(fc)] 2 (13) 

^mode 

where iV mo( j c denotes the total number of independent 
Fourier modes which fall into the (fc, /x) bin over the 15 
realizations. Note that in estimating the best-fit a v for 
dark matter, we simply set b(k) = 1 and P c = 0. 
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TABLE II: Summary of the halo catalogs. The minimum, maximum and mean mass (M m i n , M ma x and Mh) are in units of 
h~ 1 MQ, while the halo number density (nn) is in /i 3 Mpc -3 . The bias parameter, bo, is defined in Eq. (|15[1 . See Sec. IIV Al for 
more detail. 



Sample 


bin 1 (light) 


bin 2 


bin 3 


bin 4 


bin 5 


bin 6 


bin 7 


bin 8 


bin 9 (heavy) 




1.77 x 10 12 


2.49 x 10 12 


3.54 x 10 12 


4.98 x 10 12 


7.09 x 10 12 


1.00 x 10 13 


1.42 x 10 13 


2.01 x 10 13 


2.84 x 10 13 


M max 


5.54 x 10 12 


1.02 x 10 13 


1.74 x 10 13 


2.66 x 10 13 


4.04 x 10 13 


6.76 x 10 13 


1.19 x 10 14 


2.08 x 10 14 




M h 


2.96 x 10 12 


4.65 x 10 12 


7.08 x 10 12 


9.37 x 10 12 


1.47 x 10 13 


2.18 x 10 13 


3.21 x 10 13 


4.63 x 10 13 


7.03 x 10 13 


rib. 


1.57 x 10~ 3 


1.26 x 10~ 3 


9.46 x 10~ 4 


6.87 x 10~ 4 


4.87 x 10~ 4 


3.47 x 10~ 4 


2.43 x 10~ 4 


1.64 x 10~ 4 


1.09 x 10~ 4 


bo 


1.08 


1.16 


1.25 


1.35 


1.47 


1.62 


1.80 


1.99 


2.26 



Finally, in discussing the goodness of fit between the 
streaming and TNS models [i.e., Eq.([7]) with and without 
corrections] , it is convenient to introduce the reduced chi- 
squared statistic: 



Xrcd 



[Ps -body i f^j ) -fmode 



(14) 



where the power spectra -Fk-body an d -Pmodci are respec- 
tively obtained from the N-body simulations and the an- 
alytic models with the best- fit value of <r v . The quantity 
v is the number of degrees of freedom, and depending on 
the range of fitting, wet set v = 204, 364, 524, and 684 
for the maximum wavenumber used to fit, fc max = 0.08, 
0.12, 0.16, and 0.2/iMpc -1 , respectively. 



IV. RESULTS 

In this section, the results of the power spectrum mea- 
surement are presented, and a detailed comparison be- 
tween N-body simulations with analytic models is made. 
We first address the real-space clustering for dark matter 
and halos, and measure the halo bias in Sec. IIV Al We 
then move to the discussion on the redshift-space power 
spectrum, and the clustering anisotropics caused by red- 
shift distortions are shown in two-dimensional plane in 
Sec. IIVBI In Sec. IIV C| applying the multipole expan- 
sion to the anisotropic power spectra, the lowest three 
multipole spectra, i.e., monopole, quadrupole, and hex- 
adecapole spectra, are quantified and compared with an- 
alytic models. While we mainly analyze the dark matter, 
and light and heavy halo catalogs, we also examine the 
other halo catalogs in Sec. IIVDI and study how well the 
analytic models can reproduce the N-body results when 
we vary the halo mass and/or fitting range. 



A. Real-space clustering 

Since the validity and precision of analytic models 
given in Eq. (O heavily rely on the accuracy of the pre- 
diction in real space, it is important to first check the 
PT treatment, and to compare the PT prediction with 
N-body simulations. Top panel of Fig. [1] shows the ratio 



of the matter power spectrum to the smooth reference 
spectrum, P m (k)/ P nw (k), where the reference spectrum 
Pnw(k) is calculated from the linear transfer function of 
no- wiggle approximation by Ref. [53]. While the dashed 
line represents the linear theory prediction, the solid line 
is the result based on the improved PT calculation by 
Refs. [H, HI], including the next-to- leading order non- 
linear corrections up to the second-order Born approxi- 
mation. At z = 0.35, the development of non-linear grav- 
itational clustering is bit significant, and even the BAO 
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closure approx. (2nd Born) 



-+-H- 



I I I I 

bin 9 (heavy) 




Eq.(15) 



bin 1 (light) 



-+-H- 



I i I i | i I i I i I i I i | i I i I i I i 
bin 1 (light) 

° 8 % S 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 g 8 | 



bin 9 (heavy) 
i.i 



0.00 0.05 0.10 0.15 

k [h Mpc" 1 ] 



0.20 



FIG. 1: Ratios of the real-space power spectra for dark matter 
and halos. Top: matter power spectrum divided by the no- 
wiggle linear power spectrum in [53]. The symbols are N-body 
data, while the lines are analytical predictions (solid: Closure 
Approximation up to the 2nd order of the Born approximation 
[3a ]. dashed: linear theory). Middle: bias parameter defined 
by the ratio of the halo-matter cross power spectrum to the 
matter auto power spectrum [symbols: N-body data, lines: 
fitted results of Eq. (1151) ]. The results for the halo samples of 
bin 1, 2, and 9 are shown from bottom to top. Bottom: 
residual noise power spectrum, P € (k), divided by the Poisson 
noise (bin 1, 2, and bin 9 from top to bottom). 
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scales at k < 0.2/iMpc _1 exhibit a prominent enhance- 
ment in the power spectrum amplitude. Nevertheless, 
the improved PT prediction reproduces the N-body re- 
sult fairly well, and the agreement between simulations 
and PT prediction is mostly within the errorbars esti- 
mated from the 15 realizations of N-body simulations. 
This indicates that the improved PT by Refs. (3f| [Hf 
can be reliably applied to the modeling of redshift-space 
power spectrum at k < 0.2/iMpc" 1 , where the acoustic 
signature is still visible. 

In addition to the PT treatment, the halo clustering 
bias measured in real space is also an important building 
block to precisely model the redshift-space power spec- 
trum of halos. Middle and bottom panels of Fig. [1] quan- 
tify the halo bias properties based on the expressions (|T0|) 
and (ITT1) . Middle panel plots the halo bias parameter 
b(k), while bottom panel shows the residual noise spec- 
trum divided by the Poisson noise Pp i s = where 
rih is the mean number density of halos in each subsample 
(see Table |n|. The halo mass dependence is clearly seen 
not only in the halo bias but also in the residual noise 
spectrum. As increasing the halo mass, while the bias 
parameter increases, the residual noise relative to the 
Poisson error monotonically decreases and manifests a 
non-Poissonian feature, i.e., P e /Pp is 7^ 1- Although the 
mass dependence of the bias parameter is qualitatively 
and even quantitatively well-understood based on the 
halo model prescription (see 54] for a review), the non- 
Poissonian feature of the residual noise indicates that the 
origin of the residual noise might not be simply explained 
by the halo sampling process. Recently, an attempt to re- 
produce the non-Poissonian feature of the residual noise 
has been made, and the mass dependence of the residual 
noise is shown to be explained by the halo model [55$ . 
Though we do not discuss at all a quantitative aspect of 
this, the non-Poissonian feature of the residual noise po- 
tentially affects the power spectrum estimation, and the 
understanding of it is practically important to extract a 
pure clustering signal. 

Apart from the origin of the residual noise, it is worth 
noting that the scale-dependence of the halo bias as well 
as the noise power spectrum is quite moderate over the 
scales of our interest. This partly validates our treatment 
in computing the power spectrum from the TNS model 
(O (see Sec. [TTJ) . The weak scale-dependence of the halo 
bias can be well-described by a simple fitting function 
used in the literature 



b\k) 



1 + Qfc 2 
1 + Ak 



(15) 



As a reference, the fitted results of Eq. (|15p are shown 
in solid lines in Fig. [TJ and the best-fit values of the pa- 
rameter bo, which represents the clustering bias in the 
large-scale limit, are listed in Table [TT1 [67]. The resultant 
best-fit value b of the heavy sample is quite close to the 
observed value of the luminous red galaxies in the SDSS 
DR7. 



B. Redshift-space clustering in two dimensions 

We are in position to discuss the clustering proper- 
ties of dark matter and halos in redshift space, and to 
examine the validity of the analytic model prescription. 
In this subsection, we particularly focus on the redshift- 
space power spectrum in two dimensions, and compare 
the model predictions with N-body simulations. 

Let us first show Fig. [H in which we plot the ratio of 
the power spectra defined by 



P{k,n) 



{b(k) + f^} 2 P aw (k) 



(16) 



for the dark matter (top), light (middle), and heavy 
(bottom) halos. The results are shown in the (k±,kn) 
plane, where k± and fen are the wavenumber perpendicu- 
lar and parallel to the line-of-sight direction, and are re- 
lated to the quantities k and jj, through k± = k(l — /x 2 ) 1 / 2 
and fcii = k /z. Note again that the quantities P nw and / 
mean the smooth reference linear spectrum based on the 
no- wiggle approximation by Ref. [53| . and growth-rate 
parameter defined by / = dhiD + /dhia, respectively. 

In each panel of Fig. [21 the results obtained from the 
N-body simulations, the analytic models (J7J with and 
without corrections (indicated by 'TNS' and 'streaming', 
respectively) are presented from left to right. Here, the 
analytic model predictions are plotted adopting the best- 
fit values of <r v , which are derived based on the Lorentzian 
form of the Finger-of-God damping. In all of the pan- 
els, a ring-like structure originated from the BAOs is 
clearly manifest, and the acoustic feature is rather promi- 
nent along the fcj_ axis, where redshift distortions have 
no effect. On the other hand, in the line-of-sight direc- 
tion, the appreciable reduction of the amplitude is found 
in the dark matter power spectrum, while no such ef- 
fect is apparent in the halo power spectra. Rather, the 
heavy halo subsample exhibits a bit large enhancement 
at fcii > 0.1/iMpc -1 , and the effect becomes significant 
as increasing the wavenumber. Note that with the ratio 
defined by Eq. (fT6|) . the linear Kaiser effect, which also 
enhances the clustering amplitude, is effectively elimi- 
nated in Fig. [H While a damping seen in the dark mat- 
ter clustering would be due to the Finger-of-God effect, 
the absence of damping or a large enhancement in the 
halo clustering indicates that the Finger-of-God damp- 
ing seems ineffective, and/or the new effect of redshift 
distortions additionally arises and it may compensate for 
the Finger-of-God damping. Since the clustering ampli- 
tude is rather sensitive to the halo subsamples, this new 
effect seems to have a strong dependence on the cluster- 
ing bias. 

Now, looking at the analytic model predictions 
with and without the corrections, there appears a no- 
ticeable difference between the two models in the result 
of heavy halo subsample. Remarkably, the TNS model 
[i.e., Eq. (J7J including the corrections] reproduces a large 
enhancement seen in the N-body simulations quite well. 
Since the main distinction between the two models comes 
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FIG. 2: Power spectra in redshift space normalized by [b(k) + ffj, 2 ] 2 P uv ,(k) in (kx,kn) plane. We adopt the function b(k) in 
the normalization factor for halos directly measured in real space by Eq. (|10[) . while we set b(k) = 1 for dark matter. Results 
for dark matter, light and heavy halo subsamples are shown from top to bottom panels: measurements from the N-body 
simulations (left), best-fit TNS model (middle), best-fit streaming model (right). Note that in computing model predictions, 
we adopt the Lorentzian form of the Finger-of-God factor for the model predictions, and the velocity dispersion a v is determined 
by fitting the model predictions to the corresponding N-body data shown in the left panels. 



from the presence of additional terms given in Eqs. (|4]) 
and ([5]), the enhancement found in the N-body simula- 
tions can be interpreted as a non-linear effect of redshift 
distortions, characterized by the correction terms. In- 
deed, the corrections originated from the A and B terms 
exhibit a strong dependence on the halo bias parameter 6, 
and this fact naturally explains the sensitive dependence 
of the clustering enhancement on the halo subsamples. 
These findings are basically consistent with the recent 



claims by Ref. [H^, [H3] (see also Ref. (Hf), and the present 
paper further reveals an impact of the non-linear redshift 
distortions on the power spectrum in two dimensions. 

In Fig. |3l for a more quantitative comparison of the 
analytic models with N-body simulations, we fix the 
wavenumber k, and plot the ratio (|16l) as a function 
of directional cosine /z; dark matter, light and heavy 
halo subsample, from top to bottom. In each panel, the 
N-body data are shown in symbols, and the results of 
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best-fit analytic models with and without corrections, 
indicated by 'TNS' and 'streaming', are depicted as the 
solid and dashed lines. For reference, the fitted results 
of velocity dispersion er v and the goodness of fit x 2 e d are 
also shown. Compared to the dark matter clustering, a 
large enhancement seen in the heavy subsample is sta- 
tistically significant, and the discrepancy between simu- 
lations and the streaming model is clearly evident. Even 
in the absence of Finger-of-God damping (i.e., cr v =0), 
the streaming model cannot reproduce the N-body trend. 
By contrast, the model including the corrections (i.e., the 
TNS model) successfully explains the N-body results. Al- 
though no appreciable difference of the results between 
the two models is found by eye in the dark matter and 
light subsample, the resultant goodness of fit for the 
TNS model is better than that of the streaming model, 
statistically indicating that the model including the cor- 
rections successfully describes a real physical effect of 
redshift distortions seen in the N-body simulations. This 
point will be further investigated in detail in subsequent 
subsections. 



C. Multipole expansion 

The multipole expansion is an alternative technique 
for quantifying the clustering anisotropies, and it has 
been frequently applied in the literature to extract the 
cosmological information from the redshift-space power 
spectrum [1]. Especially, the lowcr-multipole spectra arc 
shown to be powerful to measure the characteristic scales 
of BAOs (e.g., d [H HI). It is also helpful for our pur- 
pose to understand visually how well the analytic models 
can accurately reproduce the simulation results. Given 
the redshift-space power spectrum in two dimensions, 
P(k,ii), the multipole power spectrum, Pg(k), is defined 



(17) 



Pt(k) = — ^- / dpV t Qj)P{k,ri, (18) 



where Vi{n) denotes the Legendre polynomial. 

Fig.|4]shows the monopole (£ = 0), quadrupole {£ = 2), 
and hexadecapole (£ — 4) power spectra (from top to 
bottom). The results are shown for dark matter, light, 
and heavy halo subsamples (from left to right), divided 
by the smooth reference power spectrum, b 2 (k) Pi, n -w(k), 
where the power spectrum P^ nw (fc) is computed from 
Pnw(k) taking the linear Kaiser effect into account: 



with the quantity j3 defined by f3(k) — f /b(k). 

Apart from the noisy hexadecapole power spectrum, 
which is largely due to the finite sampling of the Fourier 
modes along the /i direction, both of the streaming and 
TNS models are broadly consistent with the N-body re- 
sults for the dark matter and light halo subsample. A 
closer look at the acoustic feature reveals a slight dis- 
crepancy between the simulations and streaming model 
(dashed). Further, for the heavy subsample, the stream- 
ing model fails to reproduce the enhancement of cluster- 
ing amplitude on small scales. These are all what we 
found in Figs. [5]and|31 an d are consistent with Ref. pjjjj . 
The notable point is that these discrepancies are visually 
evident even for the dark matter and light halos. That 
is, while the streaming model predicts a rather clear BAO 
signal, the actual acoustic structure seen in the N-body 
simulations seems somewhat degraded, and because of 
this, the streaming model slightly overshoots the N-body 
results at low-fc, and eventually turn to underestimate at 
high-fc. By contrast, the model including the corrections 
reproduces the N-body results fairly well, and quantita- 
tively explains a slightly smeared BAOs in the cases of 
dark matter and light halos, as well as a large enhance- 
ment in the heavy subsample. 

Carefully looking at the results of the TNS model, 
however, the oscillatory behavior seen in the quadrupole 
spectrum seems to be rather over-smeared, and the 
acoustic structure becomes featureless. Because of that, 
the visual impression for the agreement with N-body sim- 
ulations is somewhat degraded. As it has been discussed 
in Ref. [33| , this is presumably due to our heterogeneous 
treatment on the corrections A and B using the standard 
PT calculations. The standard PT is known to generi- 
cally give a strong suppression on the acoustic feature in 
BAOs. As gravitational clustering develops, it may in- 
correctly lead to a phase reversal of BAOs (e.g., [59l.l60|). 
Although the corrections A and B are basically small, at 
2 = 0.35 of our samples, the non-linearity of gravitational 
clustering is strong, and the application of standard PT 
might be subtle. Nevertheless, as we will show later, the 
goodness of fit inferred from x 2 e d f avors the model in- 
cluding the corrections, and the model broadly gives a 
good agreement with N-body simulations in all halo sub- 
samples. In this respect, the model ([7]) captures an im- 
portant aspect of redshift-space clustering, and the role 
of the non-linear coupling described by A and B terms is 
quite essential. 



D. Dependence on halo mass and maximum 
wavenumber 



(1 



|/3+i/3 2 )P nw (fc) 



Pe,nw(k) = < 



(|/3 + f/3 2 )P nw (fc) 



JE;f3 2 Pnw(k) 



£ = 

1 = 2 
£ = 4 



So far, we have focused on the specific halo subsam- 
ples, light and heavy, as well as dark matter, and ex- 
ngA plored the sample dependent properties of the redshift- 
space power spectrum. Here, we examine all the halo 
subsamples, and investigate the halo mass dependence 
on the validity of analytic models in a more quantitative 
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FIG. 3: Power spectra in redshift space normalized by [b(k) + //i 2 ] 2 P nw (fc) as a function of directional cosine ft for fixed 
wavenumbers, k — 0.1225,0.1475,0.1725 and 0.1975 /iMpc -1 (from left to right). We adopt the function b(k) in the normal- 
ization factor for halos directly measured in real space by Eq. (|10|l . while we set b(k) = 1 for dark matter. The symbols are 
measurements from N-body simulations with errorbars estimated by Eq. (|13|) , while the lines are fits by the formula of Eq. |7| 
with (solid, labeled by 'TNS(L)') and without (dashed, labeled by 'streaming(L)') correction terms. Top: dark matter. Middle: 
light halo subsample. Bottom: heavy halo subsample. Note that we adopt the Lorentzian form of the Finger-of-God factor for 
the model predictions, and the velocity dispersion <r v is determined by fitting to the corresponding N-body data. The best-fit 
values of a v as well as the resultant reduced chi-squared defined by Eq. (| 14[) for each model are shown in the panels. 
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models and N-body simulations. In each panel, varying 
the maximum wavenumber fc max for the range of the fit- 
ting, the resultant best-fit value of a v divided by the lin- 
ear theory prediction a 2 lin = J dqPu n (q) / (6n 2 ) (top) as 
well as the reduced chi-squared Xi C d (bottom) are shown 
in each analytic model. Here, in computing the analytic 
models, we examine both the Lorentzian and Gaussian 
forms of the Finger-of-God damping (indicated by filled 
and open symbols, and labeled by 'L' and 'G' in the fig- 
ure legend, respectively) . The results are all plotted as a 
function of the averaged halo mass Mh for each subsam- 
ple. 

We first notice that the best-fit value cr v and Xred are 
irrelevant for the choice of the Finger-of-God dam pin g 
function, consistent with the result found in Ref. [33| . 
This is presumably because the maximum wavenum- 
ber fcmax for the range of fitting is limited, fc max < 
0.2 ZiMpc -1 , and no appreciable effect of the small-scale 
clustering appears on the scales of our interest. On the 
other hand, the best-fit a v shows a strong model de- 
pendence. While the resultant numerical values decrease 
with the halo mass in the streaming model, the fitted re- 
sults of cr v mostly remain unchanged in the TNS model, 
and interestingly are very close to the linear theory pre- 
diction cr v , lin- Further, the best-fit value in the streaming 
model is sensitively affected by the range of fitting, and 
for fc max > 0.1 /iMpc -1 , it significantly deviates from the 
linear theory prediction, and eventually becomes zero for 
massive halo subsamples. This peculiar behavior is some- 
what counter-intuitive, and indicates the breakdown of 
the model prescription. 

Indeed, for the streaming model, the goodness of fit 
indicated by x^ed becomes worse as increasing the halo 
mass and maximum wavenumber. Apparently, even at 
the best-fit values a v — 0, the reduced chi-squared x^d 
tends to keep the similar values to those obtained in the 
TNS model, and the streaming model illusively seems to 



reproduce the N-body results as well. Note, however, 
that the goodness of fit indicated by x^ed * s mere ly a sta- 
tistical measure for the discrepancy between model and 
data, and cannot be used to judge which model is phys- 
ically plausible or not. As it has been demonstrated in 
Ref. [33[ (see also Ref. (25J), the parameter estimation 
based on the streaming model can cause a large system- 
atic bias in the best-fit values of the growth-rate parame- 
ter /, while the model including the corrections correctly 
reproduces the fiducial values fairly well within the statis- 
tical error at 1-a level. Although this demonstration has 
been made in the dark matter case, the same would be 
true for the halo subsamples. Since the velocity disper- 
sion directly measured from N-body simulations is known 
to be nearly independent of the scales and halo mass 
(e.g., JHJ EH), we conclude that the streaming model 
fails to capture some aspects of redshift-space clustering, 
and the model including the corrections can provide a 
better prescription which also gives a physically reason- 
able estimate of the halo velocity statistics. 



V. IMPLICATIONS 

The goal of accurately modeling redshift distortions is 
to provide a practically useful theoretical template, and 
with a precision power spectrum measurement, we can 
tightly constrain and/or measure the growth history of 
large-scale structure as well as the cosmological distances. 
For on-going and/or future galaxy surveys, the precision 
of the power spectrum measurement is expected to reach 
at a percent-level, and even a slightly improper model- 
ing of the redshift-space power spectrum can lead to a 
serious systematic bias in the cosmological parameter es- 
timations. Then, an important question is how seriously 
the tiny deficit seen in the streaming model can lead to 
a biased estimation of cosmological parameters. 
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FIG. 5: Best-fit parameter a v divided by the linear theory prediction cj lin = J dqPn n (q) / (Qtv 2 ) and the resultant reduced 
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Here, assuming that the model including the correc- 
tions A and B, given by Eq. (J7}, is correct, we esti- 
mate the size of systematic biases caused by the incorrect 
model of redshift distortions, particularly focusing on the 
growth-rate parameter. To do this, we use the Fisher ma- 
trix formalism presented in Ref. |25l [33[ . Provided the 
survey parameters that characterize the galaxy surveys 
such as the observed redshift z c , survey volume V, and 



number density of galaxies n g , the Fisher matrix analy- 
sis enables us to estimate not only the statistical uncer- 
tainty but also the systematic bias for the best-fit value of 
cosmological parameters around the fiducial cosmological 
model. Following the Fisher matrix calculations given 
by Ref. [25|, |33[ , we specifically consider the on-going and 
up-coming galaxy surveys, BOSS '26], SuMIRe-PFS [Hj], 
and HETDEX [271 ]. and derive the expected constraints 
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on the growth-rate parameter. Table Hill summarizes the 
properties of the three galaxy surveys used in the Fisher 
matrix calculations. 

For illustrative purpose, we assume a prior knowledge 
of the real-space power spectra, P55, P$e, and Pgg. Then, 
adopting the streaming model as a wrong prior assump- 
tion for the template spectrum, we try to estimate the 
growth-rate parameter / from the redshift-space power 
spectrum in two dimensions. Notice that in addition to 
redshift distortions, the imperfect knowledge of the late- 
time cosmic expansion leads to another appreciable effect 
of clustering anisotropics, known as the Alcock-Paczynski 
effect. Here, we marginalize over the Alcock-Paczynski 
effect characterized by the distance scales, D a{z) / D Afi_d 
and H(z)/H Rd (e.g., [IE EH, S3|), where the subscript 
Ad indicates the fiducial value. As a result, the total 
number of free parameters assumed in the parameter es- 
timation is five, i.e., /, Da(z)/ -Dfid, H(z)/ -fffid, b, and 
<7 V - For each galaxy sample obtained from the three sur- 
veys, we assume the scale-independent linear bias listed 
in Table UTT] and adopt the linear theory estimate, er v ji n , 
as a fiducial value of the velocity dispersion, er v . 

Fig.[B]summarizes the forecast result for the Fisher ma- 
trix analysis based on the power spectrum measurements 
in two dimensions. The fiducial growth-rate parameter is 
depicted as solid line, and the errorbar around each sym- 
bol represents the expected 1-a constraint around the 
biased best-fit value of /, which results from the Fisher 
matrix calculation for each redshift subsample. Note that 
the maximum wavenumber fc max used in the parameter 
estimation is chosen so as to be well within the applica- 
ble range of improved PT (see Table HTT|) . where the non- 
linear gravitational clustering is still moderate, and can 
only change the real-space power spectrum by 5 ~ 10%. 

As it is clearly seen in Fig. [51 the estimated best-fit val- 
ues of / are systematically lower than the fiducial values. 
Apparently, the deviation from the fiducial value is not 
so large at each point, and the biased estimates of the 
best-fit value are sometimes inside the 1-a error around 
the fiducial value. However, the combined result with all 
the constraints strongly disfavors the fiducial model. For 
instance, if we parametrize the growth-rate parameter by 

(e.g., aii) 

f( z ) = [n m (z)r , (20) 

we obtain 7 = 0.77 ± 0.04, which significantly deviates 
from 7 = 0.55 expected from the fiducial model. 

The reason why the theoretical template neglecting 
the corrections leads to the underestimation of growth- 
rate parameter is basically explained by the mismatch 
of the overall shape of the power spectra. As shown in 
Sec. IIV C[ on large scales, the best-fit streaming model 
slightly overtakes the amplitude of each multipole spec- 
trum in both the TNS model and N-body simulations, 
and then the model eventually turns to underestimate at 
high-fc. Note that these are the outcome of the single- 
parameter fit. In order to closely match the TNS model 
and/or N-body simulations, a simple way is to further 



adjust the parameters other than a v . Among the remain- 
ing four parameters, the bias parameter b, and distance 
scales Da and H are indeed well-constrained, and their 
contribution to the systematic bias is likely to be rather 
small. In this respect, a slight change of the growth- 
rate parameter, which controls the strength of Kaiser ef- 
fect, is the only possible way to match the result of TNS 
model and N-body simulations on large scales. Note that 
whether the streaming model overestimates or underesti- 
mates the parameter f(z) heavily depends on what scales 
one weighs in the parameter estimation. In our specific 
example in Fig.[BJ we choose rather conservative values of 
fcmax- Thus, the net effect of the systematic bias appears 
as an underestimation of /. A great emphasis is that 
even the tiny discrepancy seen at low-k (see Fig. 2J can 
lead to a serious systematic bias in the on- going and/or 
up-coming power spectrum measurements. 

Recently, the WiggleZ dark energy survey has pro- 
vided a large data set of redshift-space clustering around 
0.1 < z < 0.9 from which tight constraints on 

the growth-rate parameter have been put based on the 
streaming model. After an extensive test for various 
models of redshift distortions, the authors of Ref. [64[ 
conclude that both the streaming model combining the 
fitting formula by [65| and the TNS model provide a rea- 
sonable fit to the observed clustering data, and can be 
used to derive the constraints on the growth-rate param- 
eter. Indeed, the clustering bias of this galaxy sample 
is shown to be rather small b ~ 1, and no appreciable 
distinction between the TNS and streaming models is 
manifest within the statistical errors. In this sense, a rea- 
sonable goodness of fit in both the streaming and TNS 
models sounds rather consistent with our findings, and 
the derived results on the growth-rate parameter would 
be unbiased. However, galaxy samples with a large clus- 
tering bias b > 1 will certainly exhibit a non-negligible 
effect of non-linear redshift distortions, and a tiny deficit 
in the theoretical template can lead to a large systematic 
bias on /, as shown in Fig. |6] Hence, the suitable choice 
of the theoretical template is very crucial for future pre- 
cision measurements of the power spectrum. 



VI. SUMMARY 

We have presented the redshift-space power spectra of 
dark matter and halos measured from a large set of cos- 
mological N-body simulations, and tested the analytical 
models of redshift distortions against the N-body results. 
With a particular focus on the redshift distortion effects 
of halos, we have created the nine halo catalogs over a 
wide mass range. The resultant volume and signal-to- 
noise ratio of each catalog are roughly comparable to 
those in the SDSS DR7 LRG. 

We found that the measured halo power spectrum in 
redshift space exhibits a large enhancement in amplitude 
on large scales, and the effect becomes significant as in- 
creasing the halo mass. This enhancement cannot be 
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FIG. 6: Expected constraints on the growth-rate parame- 
ter f(z) around the best-fit values (symbols), from three on- 
going/planned surveys. The results are obtained based on the 
Fisher matrix formalism assuming the streaming model as a 
wrong prior of the redshift-space power spectrum. The fidu- 
cial cosmology is depicted as the solid line, while the dashed 
line shows the best-fit model when we parametrize as Eq. (|20[l . 



simply explained by the Kaiser effect, and the popular 
model of redshift distortions called streaming model fails 
to reproduce the N-body results. The detailed compari- 
son with N-body simulations further reveals that even for 
the less massive halos with Mh > lO 13 /i _1 M0, a small 
but non-negligible discrepancy is manifest on large scales, 
and the best-fit values of the velocity dispersion a v is 
rather sensitive to the range of fitting and halo subsam- 
ples. By contrast, the model based on the perturbation 



TABLE III: Survey parameters adopted in the Fisher matrix 
analyses. The parameters, z c , V, n g , b, and fc max , represent 
the survey volume, galaxy number density, bias and the max- 
imum wavenumber included in the analysis, respectively. The 
units are in fa _3 Gpc 3 for V, /i 3 Mpc -3 for n g and foMpc -1 for 

^max - 





Zc 


V 


n a 


b 


^max 


BOSS 


0.45 


1.1 


3 x 10" 4 


2.2 


0.15 




0.55 


1.5 


3 x 10~ 4 


2.2 


0.15 




0.65 


1.9 


3 x 10" 4 


2.2 


0.15 


SuMIRe-PFS 


0.7 


0.8 


3 x 10 -4 


1.5 


0.2 




0.9 


1.1 


3 x 10~ 4 


1.5 


0.2 




1.1 


1.4 


4 x 10" 4 


1.5 


0.2 




1.3 


1.6 


4 x 10 -4 


1.5 


0.2 




1.5 


1.7 


4 x 10" 4 


1.5 


0.2 


HETDEX 


3.0 


3.0 


2.5 x 10" 4 


2.5 


0.4 



theory description (i.e., the TNS model), which includes 
non-trivial corrections to the streaming model, gives a 
better agreement with N-body simulations for every halo 
subsample. In particular, the model quantitatively ex- 
plains a large enhancement of the power spectrum am- 
plitude in heavy subsample fairy well. These results in- 
dicates that the non-linear coupling between density and 
velocity fields induces the new effect of redshift distor- 
tions, which seems to have sensitive dependence on the 
clustering bias. The corrections in the TNS model ([7]) 
can describe the observed feature of the redshift-space 
halo clustering very well. 

To investigate how the slightly improper modeling of 
redshift distortions affects the cosmological parameter es- 
timation, we have also estimated the size of systematic 
bias based on the Fisher matrix formalism. Especially 
focusing on the growth rate parameter, we consider the 
on-going and upcoming power spectrum measurements 
from BOSS, SuMIRe-PFS and HETDEX, and found that 
even the small deficit in the theoretical template of the 
power spectrum can produce a large systematic bias in 
the growth rate parameter, and the fiducial model would 
be erroneously ruled out. This is true even if we choose 
a conservative value of fc max . Hence, the accurate de- 
scription of the redshift-space power spectrum is quite 
essential, and the model including the corrections can 
be used as a reliable theoretical template of the power 
spectrum. 

Finally, we leave the following issue as a future work. 
Throughout the analysis, we have calibrated the halo 
bias relationship directly from the N-body simulations 
in real space, and the measured result of the bias param- 
eter has been used to compute the analytic model of the 
redshift-space power spectrum, after subtracting the non- 
Poissonian residual noise. Of course, this treatment is in- 
feasible in real observations, and we must combine other 
data set in order to get the information of the cluster- 
ing bias. One plausible approach is to combine the weak 
lensing measurement. Since the weak lensing measure- 
ment directly probes the matter distribution, the cross 
correlation between weak lensing and galaxy clustering 
data can give a powerful way to simultaneously charac- 
terize the clustering bias and the residual noise contri- 
butions. Synergy with the weak lensing measurement 
is worth considering, and the methodology to extract a 
pure clustering signal should be further exploited. 
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